Universal linear-temperature dependence of static magnetic susceptibility in 

iron-pnictides 
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A universal linear-temperature dependence of the uniform magnetic susceptibility has been ob- 
served in the nonmagnetic normal state of iron-pnictides. This non-Pauli and non-Curie- Weiss-like 
paramagnetic behavior cannot be understood within a pure itinerant picture. We argue that it 
results from the existence of a wide antiferromagnetic fluctuation window in which the local spin- 
density- wave correlations exist but the global directional order has not been established yet. 

PACS numbers: 74.25.Ha, 71.27.+a, 75.30.Fv 
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The recent discovery[lj of superconductivity in 
LaFeAsOi._2;Fa; has generated strong interest on the in- 
vestigation of iron-based pnictide materials. There are 
mainly two types of materials synthesized: the rare- 
earth pnictide oxide layered systems, ReFeAsO denoted 
as "1111" and the so-called "122" systems, MFe2As2 
with M=Ca, Ba, Sr, etc. Both the 1111 and 122 par- 
ent compounds are metals and have shown a spin den- 
sity wave (SDW) ordering at T ~ 130K, accompanying 
a tetragonal-orthorhombic structure phase transition 
The fact that the parent compounds of the iron pnictides 
are antiferromagnetic (AF) has attracted lots of atten- 
tion, because of the close analogy with the cuprates. In- 
deed, in the most interesting scenario, this suggests that 
the AF correlation is intimately connected to the high Tc 
in both materials. Therefore, a deeper understanding of 
the AF correlation in the iron pnictides is of particular 
importance. The purpose of this paper is to take a first 
step in this direction. 

In order to establish a microscopic theory for these ma- 
terials, two different scenarios, starting from either the 
weak or strong coupling limit, have been proposed. The 
first one invokes an itinerant electron approach in which 
the commensurate SDW ordering as well as the structural 
transition is believed to be solely induced by the Fermi 
surface nesting 0|. In contrast, the second one em- 
phasizes an As-bridged superexchange antiferromagnetic 
interactions between the nearest and next nearest neigh- 
boring local moments of irons, which serve as the basic 
driving force for both transitions ^5^ without the critical 
involvement of the Fermi surface nesting. To distinguish 
the above two scenarios, understanding of the origin of 
the SDW ordering is the key. 

Like any ordering phenomena, one can use an order pa- 
rameter ft to describe the SDW order of iron pnictides. 
In the simplest mean-field picture, n is independent of 



space and time. Above Tsdw: \n\ = and there is no 
trace of magnetism whatsoever. At Tsdw two things oc- 
cur simultaneously: a finite \n\ develops and the direc- 
tional long range order establishes. In a more realistic 
picture n is space (and time) dependent. Above TgDw 
even though locally \n\ > 0, due to the lack of directional 
order, global antiferromagnetism is absent. In the latter 
picture the SDW transition is controlled by the onset of 
directional long range order. In the following we shall 
refer to this as " SDW moment fluctuation scenario" . 

When applying the mean-field picture to the iron pnic- 
tides, one expects normal metallic behavior with no trace 
of antiferromagnetic correlation above Tsdw- As a result 
the uniform magnetic susceptibility, Xu, should be Pauli 
paramagnetic like. The Xu for both the 1111 and 122 
compounds are shown in Fig. 1 as a function of tempera- 
ture @, 0, [1, H, ■ Interestingly they exhibit a universal 
linear temperature dependence in both the undoped and 
the F-doped LaFeAsOi_a;F^ Q compounds. Clearly, this 
is inconsistent with the mean-field approach expectation. 

In the following, we argue that this linear-T suscepti- 
bility is a strong evidence for the existence of a wide anti- 
ferromagnetic fluctuation window of local magnetic mo- 
ments. It is important to emphasize, however, that the 
metallic behavior of these compounds makes the present 
local magnetic moments not quantized as those local 
atomic moment as in a Mott insulator. Interestingly, in 
undoped or highly underdoped cuprates La2-a;Sra;Cu04, 
Xu increases linearly with temperature before reaching a 
broad peak at a temperature Tmax [ill J^^t like iron pnic- 
tides. Moreover, the experimental curves can be scaled 
onto a universal curve independent of doping. This uni- 
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versal curve agrees with the theoretical result [l 
obtained for the two-dimensional Heisenberg model with 
nearest neighbor AF coupling. 

It is important to note that there are metallic SDW 
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FIG. 1: Static magnetic susceptibility Xu vs. temperature. 
The experimental data for LaFeAsOi_2,Fa; are quoted from 
Ref.Ii), SrFe2As2 from Rei.^ as 1 and Ref.H as 2, CaFezAsa 
from Ref.0], and BaFeaAsa from Ref.Q. The arrow indicates 
the experimental date of LaFeAsOi-xF^ with x — 0.04 sam- 
ple in terms of the right scale. 



systems which also show the linear-T susceptibility above 
TsDW- The best example is chromium and some of its 
alloys[15]. In the case of pure Cr, diffusive commensu- 
rate AF magnetic scattering peak had been observed up 
to temperatures T > 2TsdWj from which a very small 
effective magnetic moment (/i = 0.16 ~ 0.28/is) can be 
exacted [l6j. This suggests that the local AF SDW cor- 
relations extend to rather high temperatures. Another 
metallic AF system that shows the above linear-T sus- 
ceptibility above Tn is Nao.5Co02 which is a poor metal 
with a Neel transition at 86K[i3|. Thus the linear-T sus- 
ceptibility clearly can not be used as evidence for quan- 
tized atomic moment as in Mott insulators. 

Put it simply, such a phenomenon just implies a non 
mean-field transition into the SDW ordered state. The 
temperature range showing linear-T susceptibility is the 
fluctuation window in the Ginzburg sense. A more ap- 
propriate way of thinking is through Ginzburg-Landau- 
Wilson theory which captures the fluctuation of the SDW 
order parameter. To mimic such a theory one can write 
down an effective lattice model of flxed magnitude spin 
moments and do statistical mechanics on it. If one takes 



a classical antiferromagnetic Heisenberg model on a non- 
frustrated two dimensional lattice, it can be shown that 
the above linear-T susceptibility exists in the tempera- 
ture range < T < Tmf with /cbTmp of order the nearest 
neighbor exchange constant [3]. 

In the present paper, we prefer to start from a quantum 
spin model and do finite temperature statistical mechan- 
ics. The following two dimensional frustrated antiferro- 
magnetic Ji — J2 Heisenberg model is assumed 
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where and denote the summations over the 

nearest and next nearest neighbors, respectively. With 
J2 > >/i/2, this model captures the (tt, 0) and (0,7r) or- 
dering tendencies of iron-pnictides. Here we assume that 
at T = OK the spins are ferromagnetic ordering along 
the x-dircction and antiferromagnetic ordering along the 
y-direction. So the lattice is bipartite and divided into A 
and B sublattices. On the A (B) sublattice, the vacuum 
state is the = S (—5) state. There are two spins in 
each unit cell. 

We then use the antiferromagnetic Dyson-Maleev 
transformation to represent the spin operators. Different 
from the variational approach used by Takahashi [l^ . 
we approximate the model Hamiltonian by keeping the 
quadratic interactions of the boson operators only. Then 
the model Hamiltonian is hermitian, and it can be ex- 
pressed after Fourier transform as |20| 
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2J2S + A, Afe 2715 cos ^ + 

2 J2'5'7fe, 7fe ~ cos kx cos and a chemical potential term 
A has been introduced to make the local magnetization 
vanish at finite temperatures. By using the Bogoliubov 
transformation, we can diagonalize the Hamiltonian as 
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where ej, — y/rff, — A^, and £g is the ground state energy 
per site. From the free energy, the chemical potential A 
is determined by the following equation 
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Moreover, when a magnetic field is applied, the above 
treatment can still be carried out by including the Zee- 
man terms, and the static uniform magnetic susceptibil- 
ity is derived as 
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length is less than a lattice constant and the moments 
become effectively free, similar to the cuprates ll|, [21 1 . 



2.0 2.4 



FIG. 2: Static uniform magnetic spin susceptibility x a 
function of temperature obtained from the Dyson-Maleev lin- 
ear spin-wave calculation for S — 1 case. 



However, compared to the static uniform susceptibility 
expression obtained from the correlation function, a fac- 
tor of i has to be multiphed.(20| 

Numerical calculations for the static uniform magnetic 
susceptibility Xu can be performed at finite temperatures. 
Surprisingly, it has been found that Xu behaves as lin- 
early temperature dependence before it reaches a broad 
peak, then it can be fit as the Curie- Weiss behavior, [ioj 
There is a clear crossover regime connecting these two dif- 
ferent regimes. Moreover, as the coupling ratio of J2/J1 
is increased, the window of the linear magnetic suscepti- 
bility becomes wider. In other words, the maximal value 
of the broad peak is also shifted as increasing the cou- 
pling ratio of J2/J1. Of course, such a treatment is just 
a qualitative description of the nonmagnetic state of this 
frustrated Heisenberg model. In Fig. 2, we present the 
numerical results of Xu at finite temperatures for S — 1 
and J2/J1 = 1.0, 1.5. For the case of J2/J1 = 1-0, the 
uniform susceptibility in the temperature range between 
~ 0.9Ji, Xu can be fit as 
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a > 0. 



(6) 



Quantitatively by taking Ji ~ J2 = 55 meV estimated 
by the local density approximation (LDA) calculation^, 
we find Xo 3 X 10~^ emu/mole, which is very close to 
the experimental values extrapolated from the linear-T 
regime in Fig. 1. 

It is noted that a wide range of temperatures showing 
linear-T susceptibility can be attributed to the Mermin- 
Wagner theorem, which says that a two-dimensional 
Heisenberg system can not order at non-zero temper- 
ature. As a result, all temperatures below the mean- 
field crossover are in the fluctuation regime. The Curie- 
Weiss-like behavior will eventually recover at higher tem- 
peratures beyond Tmax ~ •/2/fcB, where the correlation 



Of course the true iron-pnictide systems show a finite 
temperature SDW ordering transition. This is due to the 
inter-layer coupling J^. In this case, we expect the fluc- 
tuation window to lie between the mean-field crossover 
and the Tsdw- In addition, due to the presence of com- 
peting interactions (Ji and J2), it is argued that above 
the SDW ordering transition, there should be an Ising- 
like transition where the symmetry between the (tt, 0) 
and (0,7r) SDW patterns are broken. Such an tran- 
sition necessarily breaks the lattice rotation symmetry, 
and as a result can trigger the tetragonal-orthorhombic 
distortion [23, [2^ ■ The persistence of the linear-T suscep- 
tibility into the doped regime implies that the SDW cor- 
relation is strong in the superconducting samples. This 
can be used as indirect evidence for the involvement of 
antiferromagnetic correlation in Cooper pairing. 

We can not overemphasize that the above quantum 
model is merely used to mimic the Ginzburg-Landau- 
Wilson description of the SDW moment fluctuation sce- 
nario. It should not be used to implies that we believe 
quantized S = 1 atomic moment exists in the system. 

Apparently, there is coupling between the SDW mo- 
ments and conduction electrons near the Fermi surface. 
For example, the SDW transition induces abrupt changes 
of the Drude weight [1^, magneto-resistance [isj . and Hall 
coefficient 25[ . In addition, angle- resolved photoemission 
experiment has shown a change of the electronic structure 
near the Fermi energy at Tsdw[2^- These experiments 
suggest that as the SDW moment orders, magnetic scat- 
tering further gap out parts of the Fermi surface and as 
a result some itinerant carriers are lost. 

Let us now switch to the electron origin of the SDW 
order. A popular point on this issue says that the 
SDW moments form because of the Fermi surface nesting 
effect 0,11. First of all, the Fermi surface nesting is an 
"instability" concept. To be precise, in the presence of 
Fermi surface nesting even infinitely weak SDW channel 
quasiparticle scattering can open the SDW gap. For a 
strong scattering, however, nesting is not required. Ac- 
cording to the band calculations, the Fermi surface are 
not well nested by the magnetic ordering wave vectors 
(tTjO) or (0,7r). In addition, the long-ranged ordering 
moment, which is a lower bound of the preformed SDW 
moment, is about 0.365/l(b for LaFeAsO 52|, and 0.873/iB 
for BaFe2As2 [13] ■ These moments are rather big since 
they are comparable with the T = ordering moment 
0.6/iB of the spin-1/2 Heisenberg model on the square 
lattice. 

Such a large magnetization moment also rules out the 
SDW transition being the mean-field SDW moment for- 
mation temperature. If that were the case, one expects 
only the electronic states at energy fc^TsDw away from 
the Fermi energy would be perturbed. Given TgDW ^ 
130K and the band structure results, we estimate an up- 
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per bound of the ordering moment to be ^ 0.02/xb, which 
is more than one order of magnitude less than the mea- 
sured value. 

In our opinion, the SDW fluctuation moment is more 
likely due to the strong short-range repulsion between 
the electrons. For example, Ref. [5| emphasizes an As- 
bridged antiferromagnetic interactions between the near- 
est and next nearest neighbor iron electrons, which serve 
as the basic driving force for SDW moment formation 
without the critical involvement of the Fermi surface 
nesting. In addition, Ref.|2^ takes the un-nested LDA 
band structure adding moderate strong Hubbard-like 
and Hund-like interactions, and obtains a good frac- 
tion of for the SDW ordering moment in an mean- 
field theory [2^. Finally the LDA-based SDW mean- 
field calculations have yielded the ordering moment be- 
tween 2.2/xb and 2.6 ^bM- However it is typical that 
all such mean-field calculations overestimate the ordering 
moment since it does not capture the long-wave-length 
directional fluctuations. 

The experimental evidences as well as the theoretical 
considerations all lead us to conclude that the SDW mo- 
ment formation temperature for the iron-pnictides ma- 
terials should occur at much higher temperature than 
TsDW- Thus there should be a "psudogap" temperature 
for iron pnictides as well. Below such a psudogap tem- 
perature, it is appropriate to consider an effective lattice 
spin model with fixed moments such as the one given by 
Eq.(l) to describe the magnetic properties of the system. 
By comparing the energy of a variety of magnetic struc- 
tures. Ma, et. al. [s*] have estimated Ji ^ J2 to be about 
55meV for LaFeAsO and 35meV for BaFe2As2. 

In conclusion, we have argued that the universal lin- 
ear temperature dependence of the susceptibility pro- 
vides a strong evidence for the SDW fluctuation mo- 
ments with strong antiferromagnetic interactions above 
the SDW transition temperature in iron-pnictides. This 
linear susceptibility can be effectively described as the 
finite temperature behavior of a Heisenberg model with 
nearest and next nearest-neighbor AF interactions. Fur- 
ther investigations are certainly needed to put our con- 
clusion on a solid ground. 
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